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II.   NOMENCLATURE 

C„        Magnus  force  coefficient  "based  upon  freestream  velocity  and 
side  area 

cyp  acY/a£ 

C.  ,  Cp     Sideforce  correlation  constants  for  small  and  large  a 

respectively 
D         Cylinder  diameter 
L         Cylinder  length 
M         Freestream  Mach  number 

CO 

p         Dimensionless  spin  rate  (v/u) 
q         Freestream  dynamic  pressure 
r         Cylinder  radius 

IL        Reynolds  number  based  on  cylinder  length 

R         Crossflow  Reynolds  number  based  on  crossflow  velocity  and 

diameter 
U         Freestream  velocity  of  spinning  missile 
U         Freestream  velocity  of  rotating  cylinder  at  9°  ot 

V  Cylinder  surface  speed 

x         Longitudinal  distance  from  nose 

Y  Sideforce 

a  Angle  of  attack 

tu        Spin  angular  velocity 

v         Kinematic  viscosity 


III.   INTRODUCTION 

The  production  of  sideforce  (Magnus  effect)  on  spinning  bodies  at 
angle  of  attack  has  been  of  interest  to  aerodynamicists  since  the  l800's 
when  G.  Magnus  published  the  first  experimental  results  of  the  drift  of 
musket  shot  due  to  spin.  Since  that  time,  it  has  been  recognized  that 
accurate  prediction  of  Magnus  forces  was  necessary  to  calculate  the  tra- 
jectory and  stability  of  spinning  shell  and  missiles.   In  spite  of  all  the 
time  and  effort  expended  since  the  l800's,  theoretical  predictions  of 
Magnus  force  production  are  extremely  crude.  Although  the  basic  mechanism 
for  Magnus  force  production  is  recognized,  the  mathematical  model  required 
for  solutions  considering  bodies  with  even  reasonably  simple  shapes  is  too 
complex  and  time  consuming  even  for  today's  sophisticated  fluid  mechanics 
calculation  methods.  The  engineer  is  forced,  then,  to  use  the  very  simple 
phenomenological  models  developed  to  date  in  conjunction  with  experimental 
correlations  to  make  estimates  of  the  Magnus  effect  on  spinning  missiles. 
This  report  gives  the  details  of  such  a  study. 


IV.   MAGNUS  FORCE  PREDICTIONS 

It  is  generally  agreed  that  Magnus  force  production  is  due  entirely  to 
viscous  effects.  An  asymmetrical  boundary  layer  is  produced  on  the  missile 
surface  by  the  combined  effects  of  spin  and  angle  of  attack.  At  very  low 
angles  of  attack  the  boundary  layer  is  attached  over  most  of  the  missile 
length  so  that  sideforce  predictions  may  be  based  upon  displacement  thick- 
ness and  shearing  stress  calculations.  At  larger  angles  of  attack  the 
asymmetrical  boundary  layer  separates  causing  a  shedding  of  vorticity  into 
the  outer  inviscid  flow.  The  shed  vorticity  induces  an  equal  but  opposite 
circulation  in  the  outer  flow  producing  an  additional  circulation  sideforce 
contribution.  The  modeling  of  such  a  three-dimensional  boundary  layer 
separation  problem  is  beyond  the  capabilities  of  present  day  theoretical 
techniques  even  for  the  relatively  simple  case  of  an  incompressible, 
laminar  flow  over  a  rotating  body  of  revolution  with  a  cylindrical  after- 
body. Theoretical  predictions  of  the  effect  of  compressibility  and  a  tur- 
bulent boundary  layer  on  Magnus  force  must  necessarily  be  even  more  crude. 

At  low  angles  of  attack  and  small  nondimensional  spin  rates  the  most 
realistic  incompressible,  laminar  flow  theory  to  date  is  that  of  Kelly. 
Kelly  considered  the  outer  flow  of  an  open  ended  rotating  circular  cylin- 
der at  an  angle  of  attack.  A  perturbation  solution  of  the  boundary  layer 
equations  was  used  to  calculate  displacement  thickness  and  shearing  stress 
contributions  to  the  Magnus  force.  The  effects  of  a  forebody  and  shed 
boundary  layer  vorticity  were  not  considered.    At  very  low  angles  of 
attack  and  for  a  high  fineness  ratio  body  this  linear  theory  yeilds  a  rea- 
sonable estimate  of  the  Magnus  force. 

2  3 
For  larger  angles  of  attack  (but  still  small)  Power  and  Iverson  * 

made  an  estimate  of  the  Magnus  force  contribution  due  to  boundary  layer 

shed  vorticity.  This  theory  will  be  outlined  in  more  detail  because  the 
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sideforce  experimental  data  correlation  used  in  this  report  originated 
from  this  work. 

Following  Kelly's  lead,  the  incompressible,  laminar  flow  about  an 
open-ended  circular  cylinder  at  angle  of  attack  shown  in  figure  V-l  was 
considered.  It  was  assumed  that  the  effect  of  shed  vorticity  could  be 
modeled  by  using  a  crossflow  analogy.   The  circulation  induced  into  the 
outer  flow  was  approximated  by  assuming  that  the  crossflow  along  the  longi- 
tudinal cylinder  axis  was  similar  to  the  unsteady  solution  for  the  flow 
over  an  impulsively  rotated  cylinder  placed  perpendicular  to  the  free- 
stream  as  shown  in  Figure  V-2.  Values  of  the  time  dependent  lift  coeffi- 
cient for  the  impulsively  rotated  cylinder  where  estimated  from  the  results 

k 
of  Thoman  and  Szewczyk.   It  was  found  that  this  data  correlates  very  well 


by  the  relationship: 


C/t)  =  f  j 


U  t   /u~D* 


(IV. 1) 


where  U  is  the  freestream  velocity,  D  is  the  cylinder  diameter,  and  V  is 
the  cylinder  surface  speed.  The  crossflow  analogy  assumes  that: 


U  =  U  sin  a 

00 

t=iTcoTl?  (IV-2) 

U        sin  a 

00 

Using  the  additional  assumption  that  f (tt~)   can  be  approximated  by  a  linear 

00 

function  the  local  Magnus  sideforce  coefficient  based  on  diameter  and 
crossflow  velocity  is  given  by: 


Cv(x)  =  C  (^— ) 
Yv  '         o\sin  a/ 


X         1 
-x  tan  a/R,u  sin2 a 

1  -  e  d 


where  C  is  a  constant, 
o 


For  small  angles  of  attack  equation  (V.3)  leads  to  an  equivalent  local 
circulation  strength  given  by 


k(x)  .  IM  .  Ci  Rc3A^ 


(iv. k) 


where  C,  is  a  constant,  R  is  the  crossflow  Reynolds  number  based  on  cross - 
flow  velocity  and  the  cylinder  diameter,  and  T     is  the  surface  circulation 
strength.   The  total  Magnus  sideforce  was  calculated  by  integrating  the 
local  circulation  contributions  over  the  length  and  adding  the  contribution 
due  to  displacement  thickness  and  shearing  stress  asymmetry.  The  result 
for  small  a  and  nondimensional  spin  rate  p  =  v/U  is: 

(IV. 5) 


Cy  =  kia  +  k2o7/U 


where : 


kl  =  Kl)^J60^-  357.37(|)^ 


-  227 


•«©  i: 


J 


k„  = 


and  the  Magnus  coefficient  is  based  on  the  side  area  LD. 

Equation  (V.5)  predicts  that  for  a  given  missile  configuration,  spin 
rate,  speed,  and  altitude,  the  sideforce  coefficient  is  nonlinear  with  a. 
The  constant  C.  can  be  evaluated  by  considering  experimental  data  for 
cases  which  are  consistent  with  the  assumptions  required  for  crossflow 
analogy  application.  It  would  also  be  expected  that  C.  might  be  a  func- 
tion of  Mach  number  because  of  compressibility  effects.  Because  of  the 
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low  crossflow  velocities  at  low  angles  of  attack,  crossflow  Reynolds  num- 
bers are  low  enough  for  the  laminar  assumption  to  be  valid  even  at  high 
forward  speeds . 

Equation's  (V.5)  nonlinearity  suggests  that  the  dominant  Magnus  force 
production  mechanism  is  shed  vorticity  at  higher  angle  of  attack.   If 
this  assumption  is  made,  a  much  simpler  correlation  suggested  by  Iverson 
can  be  developed.  Returning  to  the  result  of  equation  V.3,  an  integration 
of  the  circulation  produced  local  sideforce  results  in  a  Magnus  force 
given  by: 


CY  =  Cq  p  sin  a 


1  /  -uL   , v 
1  +  r— (e    -  1) 
Lu 


where 


_  tan  a /-n  ir  .  h 
u  =  — - — /Rdu  sin1* a 


By  expanding  the  exponential  we  find  that 


C„  = 


•  3A 
tanas in  '  a 


2    2 
1  L  tan  a     1    _L|Iil  tan  a 

Id) 


3  D 


sinua  RLU 
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sin2«  RL2 


For  reasonable  Reynolds  numbers  and  high  fineness  ratios  the  Magnus  force 
may  be  approximated  by: 


V 


=  f  ih) 


5/U       Vk       -i 

tan  a  sin  '  a  R .  " 

Li 


(iv.6) 


Equation  V.6  predicts  that  all  Magnus  data  should  collapse  along  a  line  of 

slope  C  if  dCY/dp  =  C   is  plotted  against  the  correlation  parameter 

P 

\   tan  a  sin3^  afe)    R^¥. 
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V.   EXPERIMENTAL  MAGNUS  MEASUREMENTS 

Theoretical  predictions  of  Magnus  side  force  based  upon  vorticity 
shed  into  the  outer  flow  show  nonlinear  variation  with  angle  of  attack. 
Experimental  Magnus  measurements  were  used  to  compare  actual  side  force 
production  to  that  predicted  by  theory  and  to  establish  the  values  of 
the  empirical  factors  in  equations  V.5  and  V.6.  For  this  purpose  the 
model  configuration  selected  for  study  had  a  tangent-ogive  nose  followed 
by  a  constant  2.95  inch  diameter  cylindrical  afterbody.  The  afterbody 
length  was  varied  to  provide  three  missile  length-to-diameter  ratios, 
L/D,  of  9.9,  7.8  and  6.0. 

The  model  was  machined  from  aluminum  and  carefully  balanced  to  re- 
duce vibrational  loads  when  spinning.  Figure  VT-1  shows  schematically 
the  model,  its  support  system,  the  balance,  and  the  air  motor  drive 
system.  The  l/8  horsepower  Standard  Pneumatic  rotary-vane,  bi-directional 
compressed  air  motor  allowed  missle  rotational  rates  as  high  as  3000  RPM 
to  be  tested.  An  electromagnetic  tachometer  pickup  was  used  in  conjunc- 
tion with  a  Monsanto  101B  counter  to  determine  model  spin  rate. 

A  simple  beam  balance  was  used  to  measure  the  Magnus  side  force  and 
center  of  pressure  location.  Two  four-element  strain  gage  bridge  circuits 
were  mounted  on  the  beam  flexure  such  that  only  sideforce  and  Magnus 
moment  were  measured.  Air  motor  pressure  lines  were  in  the  vertical 
plane  of  the  balance  to  reduce  pressurization  effects  on  balance  output. 
A  subsequent  balance  calibration  showed  that  for  pressures  less  than  80 
psia  the  interaction  of  line  pressure  and  sideforce  measurement  was 
negligible.  Figure  VI. 2  shows  the  wiring  diagram  for  the  balance  system. 
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The  bridge  circuit  outputs  were  amplified  and  displayed  on  a  HP  7100  B 
two  channel  chart  recorder.  The  sting  balance  system  was  calibrated  by 
hanging  known  weights  at  known  distances  from  the  strain  gage  locations 
and  plotting  the  applied  moment  verses  chart  recorder  output.  The  data 
reduction  program  in  the  Appendix  was  used  to  calculate  the  aerodynamic 
coefficients.   Figure  VI-3  shows  the  model  mounted  in  the  wind  tunnel. 

Data  was  taken  at  a  constant  nominal  Reynolds  number  based  on  missile 
length  of  1  x  10  and  at  angles  of  attack  up  to  15  degrees.  Sideforce 
production  at  seven  nondimensional  spin  rates  between  +  O.k   was  measured 
at  each  angle  of  attack.  Air  motor  supply  pressure  was  adjusted  until 
the  desired  spin  rate  was  established.  Figure  VI -k   through  VI -6  show 
the  resulting  Magnus  force  coefficients  and  the  location  of  the  Magnus 
force  center  of  pressure  with  respect  to  the  missile  base.  These  figures 
show  a  nonlinear  variation  of  sideforce  coefficient  with  angle  of  attack 
at  constant  nondimensional  spin  rate.  The  Magnus  center  of  pressure  is 
seen  to  be  approximately  one  third  of  the  missile  length  from  the  base  for 
all  the  cases  tested.   Center  of  pressure  location  was  insensitive  to 
angle  of  attack  and  spin  rate.  The  L/D  =  9«9  data,  however,  show  movement 
of  the  center  of  pressure  with  both  angle  of  attack  and  spin  rate  for 
angles  of  attack  less  than  ten  degrees.  Sideforce  coefficient  at  constant 
angle  of  attack  and  spin  rate  is  dependent  upon  missile  L/D.  Figure  VI-7 
shows  this  effect  of  L/D  on  the  Magnus  sideforce  coefficient  for  a  non- 
dimensional  spin  rate  of  0.2. 

Theory  predicts  a  linear  behavior  of  sideforce  coefficient  with 
nondimensional  spin  rate  p.  Figure  VI-8  presents  the  experimental 
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data  plotted  in  such  a  way  to  test  this  result.  The  experimental  side- 
force  coefficients  were  divided  "by  the  nondimensional  spin  rate  to 

calculate  C   =  j^C  /dp.  If  sideforce  coefficient  is  linear  in  p  all  the 

P 
data  should  collapse  to  a  single  line  when  C   is  plotted  against  angle 

P 
of  attack.  Figure  VT-8  shows  that  for  a  constant  missile  L/D  the  data 

show  an  excellent  agreement  with  this  theoretical  result.  Also  shown 

on  Figure  VT-8  is  the  linear  theory  prediction  of  Kelly  .  The  data 

show  excellent  agreement  with  the  linear  theory  at  angle  of  attack  less 

than  five  degrees.  At  angles  greater  than  five  degrees  the  data  show 

nonlinear  behavior  with  angle  of  attack. 

An  attempt  to  correlate  the  data  with  equation  V-5  was  unsuccessful. 
The  data  reduction  program  calculated  the  value  of  the  correlation  constant 
CL  for  each  experimental  point.   Values  of  C  ranging  from  zero  to  three 
resulted.  The  value  of  C  increased  as  missle  L/D  and  angle  of  attack 
increased.  No  single  value  of  C  fits  the  experimental  data.  Figure 
VT-9  shows  a  comparison  of  equation  V-5  with  the  value  of  C,  taken  as 
unity.  Equation  V-5  consistently  overestimates  the  Magnus  coefficient  at 
angles  of  attack  less  than  approximately  ten  degrees  and  underestimates 
the  Magnus  force  at  higher  angles. 

Figure  VT-10  shows  the  experimental  data  plotted  against  the 
correlation  parameter  of  equation  V-6.  Excellent  agreement  with  this 
correlation  was  found.  The  data  of  the  higher  L.D  ratio  missle  shapes 
collapsed  about  a  straight  line  of  slope  6.6.  The  L/D  =6.0  data  also 
collapsed  as  predicted  but  at  a  value  of  the  correlation  constant  equal 
to  six. 
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Figure  V.8:   Magnus  Force  Coefficient  C 
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Figure  v. 9:   Magnus  Force  Comparison  with  Equation  V.5. 
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VI.  CONCLUSIONS 

A  comparison  of  the  experimental  data  in  the  previous  section 
with  theory  leads  the  author  to  the  following  conclusions: 

1.  For  the  configurations  tested  the  Magnus  coefficient 
is  linear  in  spin  rate  for  constant  R  ,  ^   and  L/D  for  nondimensional 
spin  rates  less  than  +  O.U. 

2.  Kelly's  linear  theory  can  predict  accurately  the 

variation  of  C   with  L/D  and  R_  only  at  very  low  angles  of  attack 

P 
(in  the  case  of  the  configuration  tested  the  angle  of  attack  must 

be  less  than  five  degrees.) 

3.  The  crossflow  analogy  prediction  that  Magnus  coefficient 
is  nonlinear  with  angle  of  attack  is  accurate.  The  data,  however, 
only  match   the  theory  if  the  linear  term  in  angle  of  attack  is 
neglected.   It  might,  therefore,  be  concluded  that  although  the 
Magnus  effects  due  to  the  attached  boundary  layer  growth  (displacement 
and  shearing  stress  contributions)  are  important  at  very  small  angles 
of  attack,  the  primary  mechanism  of  Magnus  force  production  is  through 
shed  vorticity  even  at  small  angles  of  attack. 

k.      The  crossflow  analogy  leads  to  an  accurate  Magnus  force 
correlation  parameter  if  L/D  is  sufficiently  high  (in  the  case  of  the 

present  data  L/D  must  be  greater  than  approximately  seven). 

5 
5.   Iverson  has  shown  that  the  Magnus  force  coefficient  for 

supersonic  freestream  Mach  numbers  also  show  good  correlation  until 

the  crossflow  is  near  sonic.  The  correlation  constant  C  is  expected 

o 

to  be  an  increasing  function  of  the  crossflow  Mach  number  until  this 
velocity  is  reached.  At  this  point  the  parameter  C  will  remain  relatively 

2k 


constant  because  of  the  inability  of  shed  vorticity  to  reach  the  outer 
supersonic  flow.  The  experimental  value  of  C  for  incompressible 
crossflow  Mach  numbers  has  been  determined  by  this  study 
to  be  6.6  if  the  coefficient  is  based  upon  side  area.  The  experimental 
data  therefore  correlate  well  with  the  prediction: 


V  -  3.3V(L/D)5/V/U 

P 


for  angles  of  attack  less  than  fifteen  degrees 
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COMPUTER   PROGRAM  FOR  DATA  REDUCTION 

READ    (5,100)      B,TCF,RK1,RK2,N 
100       FORMAT    (5E10.0) 

READ    (5,100)   X,D,PA,TF      ,R1L 
C  B    IS   DISTANCE  BETWEEN  GAGES    (IN)    ,TCF   IS   THE  TUNNEL  CALIBRATION 

C  FACTOR    (CM  H20/PSF),K1  &  K2  ARE  BALANCE   CALIBRATION  SLOPES 

C  (IN  LB/CHART   LINES ),N  IS   NUMBER   OF  DATA   POINTS 

C  R1L  IS  DISTANCE  FROM  FRONT  GAGE  TO  BASE    (IN) 

C  X  IS  MISSLE  LENGTH   (FT),D   IS  DIAMETER   (FT),    PA   IS  ATMOSPHERIC 

C  PRESSURE    (PSF),TF   IS   TEST   SECTION  TEMP.(F) 

C  CI  IS  SMALL  ALPHA  CORRELATION  CONSTANT  AND  C2   IS  FOR   LARGE  ALPHA 

DO  1     1=1, N 

READ    (5,100)   R1,R2,RPM    ,A,DP 
C  Rl  IS  FRONT  GAGE  READING    (CHART  LINES ) ,R2   IS  BACK  GAGE  READING 

C  RPM  IS  ROTATIONAL  SPEED    (REV/MIN.),A    IS   ANGLE   OF  ATTACK, DP  IS 

C  DELTA  TEST  SECTION  PRESSURE      (CM  H20) 

Y=(RK2*R2-RK1*R1)/B 

XCP=RK1*R1/Y 

XCPR=(XCP+R1L)   /(X*12.) 
C  Y   IS   SIDE  FORCE    (LB),   XCP   IS  CENTER  OF   PRESSURE   LOCATION    (IN) 

C  XCPR   IS   NONDIM  CP  LOCATION  FROM  BASE 

TR=  TF  +U60. 

RHO=  PA*  70 . 73/ ( 1716 . *TR ) 

ESB=  0.01295*  (0.868  +.*+*d/x)*d*d*x 

Q=DP* ( 1 . +2 . *ESB ) /TCF 

V=    (2.*Q/RHO)**0.5 
C  CY  IS   SIDE  FORCE  COEFF  BASED   ON  SIDE  AREA 

CY=Y/Q*D*X) 

PHAT=3 .  l*a6*D*RPM/  (60 . *V ) 

XMU=( 2 . 26967*TR**1 . 5 )/ ( (TR+198 . 72 )*10 .**8 . ) 

RL=RHO*V*X/XMU 

XLX=  X/D 

RAD=  A/57.3 

R0=  RL*SIN(RAD)   /XLX 

XK3=  XLX*(60. 97-357- 37*XLX/(RL**. 5) -227- 01*XLX*XLX/RL)/ 

l(RL**.5)*3.11+l6/^. 

SA=SIN(RAD) 
TA=TAN(RAD) 
C  CYCP  IS  THE  CORREALATION  PARAMETER  FOR   LARGE  ALPHA 

CYCP=XLX**1. 25*SA** . 75*TA/(2 .*RL** . 25 ) 
IF    (PHAT)    3,2,3 

2  CYP=0.0 
C1=0.0 
C2=0.0 
GO  TO  k 

3  CYP=CY/PHAT 

LF    (RAD)    6,5,6    • 

6  IF   (SA)   7,5,7 

7  IF  (TA)  8,5,8 
5  C1=0.0 

C2=0.0 
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GO  TO  k 
8  C1=U . *RL** . 25* ( CY/PHAT -XK3*RAD )/ ( 3 . l*a6*XLX**l . 25*RAD**1 .  75 ) 

C2=2 . *CY8RL** . 2  5/ (PHAT*SA** . 75*TA*XLX**1 . 25 ) 
k  WRITE    (6,110) 

110  FORMAT    (IH,//  19X, 'ALPHA', 16X,,PHAT',19X,,CY,,17X,,XCFR,,18X, 
l'CYP',17X,'CYCP'    ) 

WRITE    (6,111)   A,PHAT,CY,XCPR,CYP,CYCP 

111  FORMAT      (   9X,6E20.7) 
WRITE    (6,112) 

112  FORMAT    (IH,/  39X, 'L/D' ,l8X, 'RL' ,l8X, 'RC ,17X, 'CI* ,l8X, 'C2' ) 
WRITE    (6,113)     XLX,RL,RC,C1,C2 

113  FORMAT   (29X,5E20.7) 
1  CONTINUE 

STOP 
END 
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